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Quasi-projectile events from peripheral 35 A.MeV Au+Au collisions are compared 
to central Au + C, Au + Cu and Au + Au events in the same range of excitation 
energy in terms of critical partitions and critical exponents. All the different data 
sets coherently point to a value E* = 4.5 A.MeV for the apparent critical excitation 
energy. The critical exponents r, a are compatible with the observation of a liquid- 
gas phase transition for these systems. 
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I. INTRODUCTION 

Multifragmentation has been tentatively connected to a critical phenomenon [|[] since the 
first heavy ion experiments. The recent determination of a consistent set of critical exponents 
in different multifragmentation data [Q, ^ tends to confirm this hypothesis though finite size 
corrections to scaling are unknown to a large extent. 

Deviations from critical behavior have been however observed for heavy sources 
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^, 1^, 1^ formed in central collisions and the question arises whether the multifragmentation 
mechanism depends on the entrance channel. Indeed thermal scaling near a critical point 
is expected in the case of equilibrated sources, which decay independently of the way they 
were formed. High energy experiments Q observe partitions strongly governed by the 
collective energy, and the deviations of the charge distribution from a power law of exponent 
T ^ 2 can be explained by dynamical models |0. 

However, an apparent exponent r <C 2 has been observed |P, |^, § also in central collisions 
measured at intermediate incident energies where collective components, if any, represent 
a small correction to the thermal excitation energy. A tentative explanation of this phe- 
nomenon suggests that Coulomb interaction would prevent heavy charged systems to attain 
statistical equilibration |10[ . 



In the following, by comparing Quasi-Projectile events from peripheral 35 A.MeV Au+Au 
collisions and central collisions of the asymmetric Au + Cu reactions at 25 and 35 A.MeV 
incident energy, we will show that in the latter case the increase of the Coulomb strength with 
respect to a Au source leads to asymptotic partitions of early emitted fragments coexisting 
with secondary decay products. Once this sequential emission is backtraced, thermal scaling 
is observed for several moments of charge distributions in good agreement with the QP data. 

A multiparametric fit within a Fisher droplet ansatz JTI], ^] allows to analyze two addi- 
tional data samples (central Au + C and Au + Au), at excitation energies far away from the 
expected critical energy. 

Within this ansatz all the reactions consistently point to the same set of critical exponents 
and give indication that the critical excitation energy is E* = 4.5 ± 0.5 A.MeV. Only for 
the heaviest system (central Au + Au at 35 A.MeV) a deviation from thermal scaling is 
observed, but more data at lower incident energy are needed to confirm this point. 

The plan of the paper is as follows. After a short presentation in Section II of the different 
data samples collected at the NSCL National Laboratory with the MULTICS-MINIBALL 
apparatus, the scaling properties of Au Quasi-Projectiles detected in peripheral collisions of 
the 35 A.MeV Au+Au reaction will be reviewed in Section HI. The different central data sets 
are presented in Section IV, where the correlation function technique to backtrace secondary 
fission is explained, and analyzed in Section V in terms of Fisher scaling. Conclusions are 
drawn in Section VI. 
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II. EXPERIMENT 

All the measurements discussed in this paper were performed at the K1200-NSCL Cy- 
clotron of the Michigan State University. The MULTICS and MINIBALL arrays were cou- 
pled to measure charged products, with a geometric acceptance larger than 87% of Air. More 



experimental details have been reported in Ref. |13 . 

Peripheral collisions of predominantly binary character were selected for the reaction 
Au + Au at 35 A.MeV |^. A Quasi-Projectile source of nearly constant charge and excita- 
tion energy from 1 to about 8 A.MeV was reconstructed through an event shape analysis. 
The fragment angular distributions \T^\ are compatible with an isotropic emission. The 
reproduction of the experimental charge partitions by statistical model calculations |jl6| is 
consistent with a good degree of source equilibration P, [15| . 

For the reactions Au+Cu at 25 and 35 A.MeV the most central 10% of the total measured 
cross section was selected by a charged particle multiplicity cut. Well detected events (total 
detected charge larger than 90% of the total charge) were analyzed in terms of event shape 
and the selected almost spherical events {Qfiow > 60°) are hereafter presented. The same 
procedure has been followed to select central events in the Au + Au reaction, while in the 
Au + C collision at 25 A.MeV the average fragment multiplicity < M >= 2.2 (standard 
deviation 0.2) is too low to perform any shape analysis and only the total multiplicity and 
total detected charge cut have been applied. 

Mass numbers for clusters of charge 1 and 2 were measured. To determine the cluster 
mass distribution, clusters of a given charge Z were counted on an event by event basis 
and it was assumed that N{A) = N{Z). Checks were made on the influence of the mass 
estimation. Our main conclusions do not change if we estimate the cluster mass from the 
Ao-to-Zo ratio of the fragmenting system or we assume the mass of stable nuclei. 

III. QP SOURCE : SCALING OF STATIC OBSERVABLES 

Renormalization group (RG) arguments lead, near the critical point, size distributions 
scaling as |T7 



n{A,e) = qA~^f{eA") 



(1) 
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where n{A) = N{A)/Aq is the cluster distribution normahzed to the size of the fragment- 
ing system, e measures the distance from the critical point (e = {E* — E*)/E* for events 
sorted in excitation energy bins), g is a normalization constant 0, / is the scaling 
function and r, a are universal critical exponents. 




FIG. 1: Peripheral Au + Au collisions. Top panel: scaling function for the QP source with the 
ansatz eq.^. Bottom panel: mass distribution of the QP decay products in the range A < E* < 5 
A.MeV, line: power law with r = 2.12 



For peripheral Au + Au collisions the critical exponents as well as the critical excitation 
energy of the Quasi-Projectile source (QP) were extracted from the analysis of the moments 
of the charge distribution f^. The critical exponents (r = 2.12,0" = 0.64) form a consistent 
set and are close to the values of the liquid-gas universality class. Near the excitation energy 
identified as critical from the behavior of the biggest fragment @, -E* ~ 4.5 A.MeV, the size 
distribution is well described by a power law of exponent r = 2.12 as shown in the bottom 
panel of Fig.0. 

By applying eq.(|l|) to the QP events, a scaling behavior is observed (Fig.|l] top panel) 
over the whole range in excitation energy (from 1 to about 8 A.MeV) and fragment charge 
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5 < Z < 24. Since for the lightest fragments side feeding or non-equihbrium effects could 
affect the final production rate, fragments with Z < 5 were discarded. 

In recent papers |jT2|, |1^ Elliott et al. analyze the critical behavior of the size distribution 
in terms of the Fisher droplet model. In this model the vapor coexisting with a liquid in 



the mixed phase of a liquid-gas phase transition is approximated by an ideal gas of clusters. 
A scaling around the critical point, similar to eq.(|I]) is assumed, but a different form is 
suggested for the scaling function 

n{AT) = qA-^exp{^^^^) (2) 

where AG{A, T) represents the variation of the Gibbs free energy upon formation of a 
drop of A nucleons out of an homogeneous nuclear vapor at temperature T 

AG{A, T) = AAfi - coeA" (3) 

Here A/x represents the difference in chemical potential between the two phases, e mea- 
sures the distance from the critical point e = (Tc — T)/Tc, and cq is the surface energy 



coefficient. This expression can be extended to charged systems by calculating |T6[ the 
Coulomb part of the free energy 6C in the Wigner-Seitz approximation 

AGUA Z,T) = AG{A,T) - 6G = AG{A,T) - - (l - c^e^f') (4) 

where R is the radius of a fragment of mass A, Cc is a constant related to the density of 
the liquid phase [|19| and /5 is a critical exponent. 

To verify the ability of this ansatz to extract information about criticality, we performed 
a multi-parameter fit of the QP data (over 200 data points). The values of the critical 
parameters from the fit were then compared with those already known from the analysis of 
the moments of the charge distribution 

As in Ref. |12] the temperature T in eqs.(^, ^ has been replaced by ^/8E*JA. This 



ansatz has not to be interpreted as a realistic temperature estimator, since the assumption 
that close to the critical point nuclear systems behave like a degenerate Fermi gas is highly 



disputable and not supported by experimental data pO|- This Fermi gas ansatz has only 



to be considered as an effective parameterization of the scaling function which allows a 



direct and quantitative comparison with the similar analysis performed on Isis data [12 . 
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In principle A/i as well as Cq |T2[ and Cc can depend on temperature. To keep parameters 
under control we have taken the simplest parameterizations that allows a good quality fit of 
the whole cluster size distribution, cq and Cc were treated as free parameters, independent 
from the temperature, the exponent (3 was fixed at the value 1/3 as resulted from previous 
analysis on the QP source 0. Whatever is the parameterization chosen for A/x, its value at 
E* = E* is compatible with zero. 



Peripheral collisions 


T 


a 


E* {A.MeV) 


Co 




5C = 0, An fixed at 


2.05 ± 0.01 


0.66 ± 0.06 


4.19 ± 0.06 


3.3 ± 0.6 


1.5 


6C = 0, A/x = const. 


2.07 ± 0.02 


0.70 ± 0.03 


4.5 ± 0.4 


3.0 ± 0.3 


1.6 


5C = 0,Af, = f{T) 


2.10 ± 0.02 


0.66 ± 0.02 


4.5 ± 0.1 


8.0 ± 0.3 


1.3 


6C eq.(|), A^l = f{T) 


2.08 ± 0.02 


0.66 ± 0.02 


4.40 ± 0.05 


10.9 ± 0.4 


1.7 



TABLE I: Values of the critical exponents t, a, critical excitation energy E*, surface energy 
efficient Co and of the Fisher scaling fit through Eq.s^j^ 



coe 



In Table | we report the results obtained for two different parameterizations for A/i (a 
constant value or a polynomial of order two in E*/A) together with some results obtained 
switching on/off the Coulomb term of eq.^. Checks were also made on the parameterization 
chosen for the Coulomb term. Results compatible with those reported in Table | are obtained 
with the parameterization used in Ref. [T^ . 



The fit through Fisher scaling function (eq.^ and following) produces a scaling behavior 
for the data, comparable to the RG ansatz as can be seen in Fig.|^ (bottom panel). Critical 
exponents and the critical excitation energy, reported in Table |, are in agreement with the 
previous analyses of the moments of the charge distribution. It is interesting to remark 
that similarly to the RG ansatz the scaling obtained with the Fisher scaling function is 
observed in the sub-critical and in the supercritical domain, even if in the context of the 
Fisher model eq.(|]) has no physical meaning for E* > E*. This is more quantitatively shown 
in Figure |^ (top panel) which compares the measured production yields of a few clusters 



with the results of the fit. As already observed in Ref. |21] a parameterization Afi = f{T) 
allows the Fisher scaling function to reproduce the data well above the critical excitation 
energy. On the contrary, when Afi is treated as a constant, the scaling is observed, mainly 
for large cluster sizes, only up to the critical excitation energy. 



E' (A.MeV) 




FIG. 2: Top panel: Cluster yields of the QP source as a function of the excitation energy (symbols) 
and fits with eq.^. Dotted line: 5C eq.^ and A/x = f{T), solid line: 6C = and Afi = fiT), 
short dashed line: 5C = 0, A/.i = 0, long dashed line: 5C = 0, = const.. Bottom panel: 
scaled yield distribution versus the scaled temperature for the cases: 5C = 0, A^ = (left) and 
Afi = f{T) (right). To visualize the results on the same picture a constant horizontal shift C = 1 
is given to the latter distribution. 



The similarity of the resuhs obtained with the RG and with the Fisher scahng demon- 
strates that the values of the critical parameters are very robust with respect to the detailed 
shape assumed for the scaling function. On the other hand this also implies that the ad- 
equacy of the Fisher ansatz to the data is not a sufficient condition to give any detailed 
physical meaning to the value of the fit parameters Afi, cq, Cc nor to the sorting parameter 
^y8E*/A. As an example, the quality of the fit and the value of the critical parameters do 
not drastically change if one suppresses the Coulomb as well as the volume term |jl8| (see 
Table |); the value of cq strongly depends on the parameterization of the scaling function 
and the full scale of the scaled yields depends also on the level density parameter assumed 
to estimate the temperature. 
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Peripheral collisions 


T 


a 


E* (A.MeV) 


Co 




E* = 1-2 {A.MeV) 


2.5± 0.2 


0.66 ± 0.03 


4.5 ± 0.6 


10.0 ± 0.7 


1.4 


E* = 2-3 (A.MeV) 


2.1 ± 0.1 


0.65 ± 0.03 


4.5 ± 0.3 


10.0 ± 0.9 


1. 


E* = 3-4. (A.MeV) 


2.10 ± 0.08 


0.68 ± 0.05 


4.5 ± 0.2 


10.0 ± 1.5 


0.7 


E* =4-5 (A.MeV) 


2.08 ± 0.04 


0.66 ± 0.02 


4.5 ± 0.2 


10. ± 5. 


1.0 


E* = 5-Q (A.MeV) 


2.08 ± 0.02 


0.66 ± 0.06 


4.50 ± 0.02 


10.1 ± 0.2 


1.1 


E* = 6-7 (A.MeV) 


2.08 ± 0.02 


0.66 ± 0.04 


4.50 ± 0.02 


10.1 ± 0.1 


1.5 


E*>7 (A.MeV) 


2.08 ± 0.02 


0.66 ± 0.03 


4.46 ± 0.02 


10.0 ± 0.1 


1.1 



TABLE II: Values of the critical exponents t, a, critical excitation energy E*, surface energy 
coefficient cq and of the Fisher scaling fit trough Eq.s^^,^ in bins of excitation energy. 



The Fisher scahng procedure has also been performed on several intervals of the excitation 
energy and the results are summarized in Table |T|. This test allows to validate a Fisher 
scaling analysis on central collisions, where the excitation energy is distributed in a narrow 
interval. 

As can be seen in Table || for all the energy intervals the fitting procedure recognizes the 
same critical point even if the critical energy is not contained in the set of data. The critical 
exponents are again in agreement with those previously obtained. 

The uncertainty in the critical exponents and in the critical excitation energy decreases 
for increasing E*, signaling that only an approximate description of the critical region can 
be provided by events containing information only on the vapor, as previously found in the 



study performed in Ref. |^ on Lattice Gas events. Another signal is given by the value 
of Afi at E*. While this value is compatible with zero for all the excitation energy bins 
E* > 2 A.MeV, this is not the case for the first interval of excitation energy, where it results 
AiJ,(E* = E*) = 7.2 ±0.1 A.MeV. This does not seem merely related to the distance from 
the critical excitation energy : indeed "overcritical" events 2 A.MeV above E* recognize the 
critical point better than events 2 A.MeV below. 

To conclude, we have shown that the Fisher scaling technique is extremely powerful for 
the determination of the critical exponents and the critical excitation energy. However, the 
adequacy of the Fisher ansatz to the data is not sufficient to give a physical meaning to the 
other fit parameters. 
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IV. CENTRAL COLLISIONS 
A. Equilibration of the emitting sources 

We have shown in the previous Section that a critical-hke behavior is observed around 
E* — 4.5 A.MeV for QP events. To explore the effect of the entrance channel and/or the 
size of the system we will now turn to the analysis of central collisions. 

Firstly we show that the considered central events are compatible with the decay of 
equilibrated sources and therefore can be good candidates for an analysis in terms of thermal 
scaling. The relevance of the thermostatistical analysis depends indeed on the approximation 
at which an equilibrium is realized. As a general statement, the degree of approximation of 
an equilibrium is indicated by the degree of the agreement of data with statistical models 
containing the same constraints as the data. 

To check the quality of the source selection criteria, a standard procedure consists in 
verifying that events are spherically symmetric in momentum space. Of course it is very 




FIG. 3: Fragment charge Z versus the parallel velocity along the beam direction in the center of 
mass reference frame for central Au + C events at 25 A.MeV (panel a), Au + Cu events at 25 
(panel h) and 35 A.MeV (panel c) and Au + Au central collisions at 35 A.MeV (panel d). 

likely that in the dynamical preparation of the nuclear source the shape degree of freedom 
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may not be completely relaxed; such a case could in principle be addressed within a statistical 
ensemble where deformation is explicitly accounted for. However restricting the analysis to 
spherical systems guarantees that preequilibrium effects do not pollute the statistical sample. 

An observable suggested to discriminate between single source and multiple sources 
decay events, is the correlation between the velocity along the beam direction and the charge 
of the fragments. 

This observable is shown in Fig.^ for central collisions of the reactions Au + C at 25 
A.MeV incident energy (panel a), Au + Cu at 25 (b), 35 A.MeV (c) and Au + Au at 35 
A.MeV (d). In this representation one can observe that the source is at rest in the centre 
of mass reference frame and also that the general trend for fragments is to be isotropically 
emitted. 

In Figure |^ we compare the charge distributions to the predictions of statistical models. 

The Au + C data are the standard prototype of low energy compound nucleus deexcita- 
tion. Due to the strongly forward focused kinematics and to the minimum detection angle 
{Oiab > 3°), complete events are detected only in the case of fission of the compound. The 
charge distribution displayed in Figure || (top panel) is well reproduced by the Gemini ||2^ 
code, which describes fragment production as a sequence of binary fission-like emissions. 
The excitation energy needed to the model to reproduce the experimental distribution (1.4 
A.MeV) agrees with the experimental value. The small shift between the experimental and 
simulated charge distributions is presumably due to the (small) difference between the size 
of the experimental and simulated source. While Gemini can not be run for sources heavier 
than Au, in the data the charge collected by the two fission fragments is the charge of the 
completely fused system. 

The charge distribution of the central collisions extracted from the Au + Cu 25 and 35 



A.MeV reaction are compared to statistical SMM model calculations [|T6|, where the mass 
of the fragmenting source has been assumed equal to the total mass of the system and the 
freeze out volume has been taken as ^ = 3Vo. The model well reproduces the data with 
input excitation energies (3.3 and 4.1 A.MeV for the 25 and 35 A.MeV data respectively) 
close to the average energies measured by calorimetry. 

In the case of the 35 A.MeV Au + Cu sample the excitation energy is close to the value 
where the QP shows critical partitions with r ~ 2. However, as it appears from Fig.^ the 
central events distribution is much flatter and can be approximately fitted only with an 
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FIG. 4: Symbols: inclusive charge distributions for central Au+C events at 25 A.MeV, Au+Cu at 
25, 35 A.MeV and Au + Au at 35 A.MeV from top to bottom, respectively. The reported excitation 
energy and standard deviation are the experimental calorimetric values for each sample. Lines: 
filtered GEMINI (top panel) and SMM (other panels) simulations. 



effective exponent r^ff ~ 1.4. The average excitation energy is lower for the 25 A.MeV Au + 
Cu sample, however also for this reaction a two parameter power law fit of the distribution 
gives an apparent exponent Tg// ~ 1.4. 

The decay of the Au + Au source (bottom panel of Fig.^ is also very well described by the 



statistical SMM model [0, A moderate contribution of radial flow (about 1 A.MeV) 



seems to be superimposed to the thermal excitation energy and the freeze-out density of 
this system, evaluated from the backtracing of experimental data as about 1/3-1/6 of 
the normal density. The apparent exponent Tg// of the charge distribution is about 1.2 and 
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again the deviation from a power law behavior could be due in this case to the distance 



from the critical point ||26[ or to some other (size or entrance channel dependent) physical 
effects m. 




1 10 30 50 70 Z 

FIG. 5: Charge distribution of the largest fragment (triangles), the second largest (squares) and 
the remaining fragments (circles). Upper (lower) panel: central 25 (35) A.MeV Au+Cu events. 



The same kind of agreement with statistical models is also observed, for all the considered 
reactions, for the charge partitions ^ as well as for the fragment kinetic energies. The 
close similarity between statistical models and data, together with the isotropy of fragment 
emission (see Fig.|D, suggests that this set of data can be analyzed an analysis in terms of 
thermal scaling. 

To better understand the low value of r^ff for the Au + Cu central collisions, we show in 
Figure ^ the charge distribution of the biggest, second biggest and remaining fragments for 
both samples. One can see that the charge distribution of light fragments (circles) is steeper 
than the total distribution and that the flattening, resulting in Xe// ~ 1.4, is mainly caused 
by the second largest fragment, whose distribution covers a wide charge range. 

In the next Section we give an interpretation of this behavior. 
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B. Prompt or sequential emission? 



As suggested in Ref. ||2^, the observed tail towards small Z values of the second largest 



fragment could be due to secondary decay of the largest primary fragment formed at the 
freeze-out, because of the reduced fission barriers of a source heavier than the QP one. 
According to this interpretation, fragments generated later than the freeze-out time should 
keep memory of their later emission by showing a strong velocity correlation. This should 
be apparent in the two-fragment velocity correlation function, defined as 

1 + = C,p^ (5) 

J- back \JJred) 

where Vred = ("^i — Vj)/ a/ {Zi + Zj) is the reduced relative velocity of fragments i and j 
with charges Zi and Zj] Y{vred) and Yback{vred) are the coincidence and background yields for 
fragment pairs of reduced velocity Vred and C = Nf,ack / Ncoinc where Ncomc and Nback are the 
total number of coincidence and background pairs. The background yield was constructed 
by means of the mixed event technique |^ . 



This correlation function is in principle affected by a space-time ambiguity. However, if 
the density of the sources is comparable (as indicated by statistical model predictions) this 
observable measures the fragment formation time. 

Figure || displays the correlation functions of the reduced velocity for all the sets of central 
events. These observables are also compared to the correlations obtained from QP events 
at the same excitation energy as the central sets. 

For central Au + C events at 25 A.MeV the correlation function is very similar to the QP 
one, especially when only fission events are considered (dotted line). The peak at Vred ~ 10 
(c/1000) corresponds to the Coulomb reduced velocity of two touching charged spheres. 

For central Au + Cu events the correlation function of the reduced velocity between the 
two heaviest fragments in each event signals emission times longer than in the QP decay at 
the same excitation energy. The Coulomb hole is narrower than in the QP case and the 
peak corresponding to the reduced velocity of two touching charged spheres, observed in the 
Au + C reaction, is still present even if the excitation energy of the source is larger. 

For the Au + Au central collisions the correlation function of the reduced velocity (Fig.|^ 
bottom panel) of the two heaviest fragments does not give indications of a late emission. 

Several interpretation of this behaviour can be given. One is that for Au + Au central 
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FIG. 6: Correlation functions of the reduced velocity for the two largest fragments in each event. 
From top to bottom: central Au + C events at 25 A. Me V, Au + Cu at 25, 35 A.MeV and Au + Au 
at 35 A.MeV. Triangles: correlation between detected fragments, open symbols: the same correla- 
tion after re- clusterization for central events(see text). Dashed lines: QP correlation at the same 
excitation energy as the corresponding central events. 



events all the detected fragments were originated at the freeze-out time. Another alternative 
interpretation is that the increased Coulomb strength with respect to a Au source also affects 
the decay of fragments other than the largest one. 

More sophisticated correlation techniques would be needed to answer this question. 

For the Au + Cu reactions we can use the correlation function information to backtrace 
the sequential emission. We reclusterize the two heaviest fragments with a Gaussian shaped 
probability distribution, centered at the Coulomb reduced velocity and with a width given by 
the left tail of the correlation function. When the two heaviest fragments are reclusterized, 
they are replaced in the event by a primary fragment with a charge equal to the sum of their 
charges and a velocity vector given by the weighted mean of their velocities. As shown in 
Figure after this procedure the correlation function has a very similar shape to the QP 
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one in the same excitation energy range. 

C. Critical exponents 

After the backtracing of secondary fission the cluster size distribution for the central 
collisions Au + Cu at 25 and 35 A.MeV are in very good agreement with the QP size 
distributions at the same excitation energy (see Fig.|^. In particular the 35 A.MeV reaction 




20 40 60 80 

A 

FIG. 7: Cluster size distribution for central Au + Cu events (open symbols), after re- clusterization 
of the two heaviest fragments with the procedure described in the text. The line represents the power- 
law from Eq.s^^,^ describing the QP events at the same excitation energy as central Au + Cu 
25 A.MeV collisions (top panel) and Au + Cu 35 A.MeV collisions (bottom panel). 

follows a power-law with t ^ 2 a.s appears from the bottom panel of Fig.^. The fact that 
at similar excitation energies the central and peripheral sources show a similar fragment 
distribution strongly supports the thermal scaling. 

In addition this means that, at least up to a source charge around 100, the increase of 
the Coulomb in the entrance channel only affects the decay of primary fragments, but not 
the bulk multifragmentation. 
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A power-law with exponent r fa 2 corresponds to partitions with very large size fluctua- 
tions [0, as signaled by the Campi plot of the central events (Fig.|D, showing a "liquid" as 
well a "gas" branch. 

1 







-3 -2 -1 

lln«S2>/Z/"^/') 

FIG. 8: Campi plot for the QP source (dashed line) and for Au + Cu 35 A.MeV central collisions 
(symbols). In the left panel the unsealed distribution is given, in the right panel scaled moments 
are shown. 



It has to be noted, however, that a bimodal shape of the Campi plot is a necessary, but 
not sufficient, signal to claim the observation of a phase transition. Indeed simulations. 



which do not contain criticality, but only take into account the phase space |18] can produce 



Campi-plots apparently similar in shape. Only critical exponents allow to identify 
criticality. 

A strong check of compatibility among QP critical events and Au + Cu central events 
consists in the finite size scaling of the moments of the size distribution. 

In Ref. it has been shown that, as for percolation events, in finite nuclear systems 
near the critical excitation energy the size of the largest cluster and the second moment 
scale for the size of the system following the hyper-scaling relation [|T^, . 

Following Ref. [^] the largest cluster size is expected to scale as S^^^, where S is the 
size of the system, D = is the fractal dimension, ci = 3 is the Euclidean dimension. The 
second moment 5*2 is expected to scale as S^'^"'^^. 

In the right panel of Fig.^ the correlation between the size of the largest fragment and the 
average second moment is shown for the central source, scaled ||32| via the critical exponents 
of the QP source [Q. The observed agreement shows again the validity of the scaling |33 
and the compatibility of critical exponents for the two systems. 
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In this section we have presented "compatibihty" checks between the static observables 
of emitting sources, different in size and in the way they were formed. However, our aim is 
to extract a complete set of critical exponents also for sources formed in central collisions. 

The procedure adopted for the QP in Ref. cannot be applied to central collisions, 
since the excitation energy is distributed in a very narrow interval. However we have shown 
in Section HI that even for sets of data that do not contain the critical point it is possible 
to get information on critical excitation energy and critical exponents by working with the 
Fisher scaling function eq.(||). 



V. FISHER SCALING FOR CENTRAL COLLISIONS 

In order to get critical parameters for the central events, independently of the information 
extracted from the QP data, the fitting procedure described in Sect. HI has been applied to 
the four sets of central events discussed above. 



Central collisions 


r 


a 


El {A.MeV) 


Co 




Au+C 25 A.MeV 


2.0 ± 0.6 


0.64 ± 0.02 


4.5 ± 0.1 


10.0 ± 0.2 


7.8 


Au+Cu 25 A.MeV 


2.10 ± 0.04 


0.60 ± 0.02 


4.4 ± 0.1 


10.0 ± 0.6 


1.8 


Au+Cu 35 A.MeV 


2.05 ± 0.04 


0.68 ± 0.03 


4.5 ± 0.2 


9. ± 1. 


1.3 


Au+Au 35 A.MeV 


2.1 ± 0.1 


0.68 ± 0.03 


4.51 ± 0.01 


10.0 ± 0.1 


1.7 


All central 


2.05 ± 0.02 


0.66 ± 0.02 


4.40 ± 0.03 


13.0 ± 0.4 


2.0 



TABLE III: Values of the critical exponents t, a, critical excitation energy E^, surface energy 
coefficient cq and of the Fisher scaling fit trough Eq.s^,^,^. 



The resulting scaling functions are shown in Figure ^ and the corresponding critical 



parameters are reported in Table |H . 

Also the data samples not covering the critical region point to the same critical energy 
and to compatible critical exponents. The fit is of high quality and the values of Afi at E* 
are compatible with zero, within the error, except for the compound nucleus de-excitation 
events {Au + C reaction). In this case the value of AyU at is large (2.9 ± 0.06 A.MeV) 
and the quality of the fit is poor. For this reaction, however, the excitation energy is very 
far from the critical one. As found for the QP events (see Section HI) and in the study 
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-1 -0.5 0.5 1 

A" £/T + C 

FIG. 9: Fit of central events with the scaling function ansatz eq.lj^. To visualize the results on 
the same picture a constant horizontal shift C is given to the distributions. 

performed on a Lattice Gas |21|, an approximate description of the critical region can be 
provided by events containing only information on the vapor. 

To check whether the scaling function is able to describe the decay of finite charged 
nuclear matter, independently on the size of the source, a subset of events of equal statistics 
has been selected for the four sets of central events (1000 events for each reaction) and Fisher 
scaling (H) has been repeated on this sample. The resulting scaling function is displayed 
in Fig.0 and the parameters are again reported in Table |ITl| (label : 'AH central'). The 
good quality of the scaling and the perfect agreement of the critical parameters with the 
Quasi-Projectile results are a striking proof of the universality of the scaling: the critical 
behavior and the whole shape of the scaling function are independent both on the size of 
the system and on the entrance channel. This implies that the Fisher scaling technique is 
extremely powerful and the determination of the critical parameters is robust with respect 
to the possible effects induced by the Coulomb interaction and by the reaction mechanism. 

A word of caution is however necessary. In Section II we have presented a comparable 
quality scaling obtained with the RG ansatz (Q) which applies to generic critical phenomena 
like percolation, not necessarily connected to a thermal transition. It is also well known that 
the liquid-gas critical exponents are very close to the percolation ones. Therefore other non 
static observables have to be analyzed before one can identify the transition as belonging to 
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the liquid-gas universality class, but this is beyond the aim of this paper ^6 . 
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10 



-5 



-0.4 



-0.2 







0.2 



0.4 



0.6 



0.8 



A"£/T 



FIG. 10: Simultaneous fit of an equal statistics subsample of the central events with the scaling 
function ansatz eq.^). Symbols as in Fig^. 



The decay of the Quasi-Projectile formed in peripheral 35 A.MeV Au + Au collisions has 
been compared to different sets of central events in terms of critical partitions and critical 
exponents. 

All reactions are compatible with a critical behavior at an excitation energy of about 4.5 
A.MeV, with critical exponents for the size distribution t ^ 2,a ^ 0.65. Analyzing the 
velocity correlations between the two heaviest fragments we have demonstrated that fission 
is an important deexcitation channel for the secondary decay of heavily charged systems 
{Z source ~ 100). If this dccay is back-traced we observe a perfect scaling with the size of 
the source between the peripheral QP events {Zsource ~ 80) and the central Au + Cu ones 
{Z source ~ 100). On the other hand the size distribution of central Au + Au events is flatter 
than the QP one at a compatible deposited energy, and this fact cannot be ascribed to 
secondary fission. 

The impressive quality of the scaling and the compatibility of the results of central and 
peripheral collisions does not demonstrate unambiguously the thermal nature of the tran- 
sition, though it certainly is a strong circumstantial evidence for its existence. Indeed the 
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shape of the scahng function (i.e. the parameters Cq and A/i) depends non neghgibly on the 
technical assumptions made in the procedure of the fit and a similar scaling is obtained with 
the ansatz Eq.([l|) that does not assume a thermal nature for the transition. This means 
that the different size distribution (implying a different behavior for Afi) observed for the 
Au + Au system can be ascribed to a different thermodynamical path for the state variables 
explored in central collisions of symmetric systems, but also to non thermal effects due either 
to the entrance channel (signaled by the onset of collective flow @]) or to the high charge of 
the system (Coulomb explosion for too heavily charged sources [0, §)• 

The scaling "per se" does not determine the order of the transition either. Indeed in 



many different statistical calculations ^ |35| first order phase transitions in finite 



systems lead to observations typical of continuous transitions including critical behaviors 



and critical exponents. Further work to assess this point is in progress [36| 
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